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ABSTRACT 

A quantitative theory of ionization in diffuse clouds is developed 

which includes H + charge exchange with 0. Dissociative charge exchange 

+ + + 

of He with plays an important role in the densities of H and He . 

The abundance of HD is also discussed. 
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I. Introduction. 

The densities of electrons and ions play an important role in determining 
the physical and chemical properties of interstellar clouds. In this article 
we present a theory for the ionization of interstellar clouds which should he 

21 -2 

applicable to diffuse and moderately thick clouds (column densities N £ 7 x 10 cm ). 
This work is an extension of our recent discussion of ionization in diffuse 
clouds (Glassgold and Langer 1974* especially Appendix A). 

There are two main sources of electrons in interstellar clouds: the ioniza- 

tion by the radiation field of heavy trace elements with ionization potentials 
I < 13.6 eV (e.g. carbon), and cosmic ray ionization of atomic and molecular 
hydrogen and helium. In addition to radiative recombination, ions are removed 
by charge exchange, dissociative recombination, and ion-molecule reactions. 

The last two types of reactions become important as soon as the interstellar 
clouds contain an appreciable amount of molecular hydrogen. Charge exchange 
and ion-molecule reactions only affect the electron density indirectly, in that 
the transformation of one ior. into another affords the possibility of a 
different recombination process. In particular, the rate constant for 

1 

dissociative recombination is many orders of magnitude larger than that for 
radiative recombination. Field and Steigman (1971) pointed out that charge 
exchange of H + with atomic oxygen is important in determining the ionization 
of the interstellar medium. When molecular hydrogen is present in interstellar 

* 4 * + 

clouds, the 0 produced by charge exchange interacts with Hg to form OH , and 
in subsequent interactions with and electrons a variety of oxygen bearing 
molecules and molecular ions are formed. The primary effect of these reactions 
is to remove H + and electrons at a greater rate than if only radiative recombina- 
tion were operative (Dalgarno, Oppenheimer, and Berry, 1973; Watson, 1973; and 
Black and Dalgarno, 1973). The exact rate depends sensitively on temperature 
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since 0 + H + charge exchange is endothermic by 232°K. Some of the molecules 
formed from the chain beginning with 0 , namely OH and H 2 0, can c ^ ar & e exchange 
exothermically with H and He , but their effect can be neglected if n(0H)/n(0) 
and n(H 2 0)/n(0) < lCf 2 . 

We also include a reaction which has not been previously considered in 
problems of this type, dissociative charge exchange of He and H^, 

He + H 2 H + H + He ( 1 ) 

The rate constant for this reaction has recently been measured as 

* Kf^cm^s"" 1 (Biondi, 1974). Despite this low value, reaction (1) is 

+ 

an important source of H when the fractional abundance of H 2 is high 
(f > 0.6) and it can dominate the recombination of He even for small 
values of f . The parameter f is defined in terms of the atomic and molecular 
hydrogen densities, 

f * 2n(H 2 )/n n * n(H) + 2n(H 2 ) . (2) 

t + 

Deuterium also can charge exchange with H , 

D + H + -v H + D + , (3) 

but this does not lead to a significant reduction in n(H + ) since D + reacts 

+ 

rapidly with H and H 2 to produce H again, 

D + + H + H + + D (4) 

D + + H, -*• HD + H + .... , (5) 

The HD density, n(HD), can be used to estimate the primary cosmic ray 
ionization rate on hydrogen, £ , through its connection to H* via 
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reactions (3) and (5) (Black and Dalgarno, 1973; Jura, 1974; 0 1 Donnell and 
Vat son, 1974). 

In this article we will discuss quantitatively the role charge exchange 

+ 

and related ion-molecule reactions play in the electron, E , and HD abundances 
for diffuse and moderately thick clouds. The results will depend on a number 

of physical parameters of the cloud, such as temperature T, density n, molecular 

+ m 9 

hydrogen fraction f, and heavy element abundances £. The H concentration is 
crucial for the production of interstellar molecules through gas-phase ion- 
molecule reactions. The electrons are important in determining the recombination 
of all ions and thus the abundances of all species which are substantially 
ionized. Finally, the electron density is also basic in the thermal balance 
of the cloud because electrons colli sionally excite and de-excite cooling 
transitions . 

II. Calculation 

The electron and ion abundances will be determined by solving the 
chemical balance equations for all the constituents, using the constraints 
of charge neutrality and conservation of number density for each nuclear 
species. The problem is very complex, hut it can be solved approximately 
in a number of different limits. We restrict ourselves to a regime 
characteristic of many of the clouds which have recently come under 
study, namely diffuse to moderately thick clouds with E(B - V) £ 2. 

These clouds are relatively cool (T £ 100° K) and contain substantial 
amounts of molecular hydrogen (f > 0.1 ). Extensive observations have 
recently been made on diffuse clouds with the Princeton UV spectrometer 
on the Copernicus spacecraft. The column densities deduced indicate 
that f ■ 0.1 is reached when the total hydrogen column density into the 



cloud, N = N £h]) + 2N £h 2 J , becomes ^ 5x lO^cnT 2 . This result is supported 
by theoretical studies (Hollenbach, Werner, and Salpeter, 1971; Glassgold and 

r 

Langer, 1973a and 1974; Jura, 1974). Kinetic temperatures in the range from 
70 - 120°K have been determined for these clouds from the relative populations 
of the lowest J = 0 and 1 states of (Spitzer and Cochran, 1973; Spitzer, 
Cochran, and Hirshfeld, 1974). The temperatures in somewhat thicker clouds 
may be expected to be lower because of the attenuation of external heating 
radiations. Thus we are interested in determining the electron and ion 
abundances in regions where the interstellar gas is predominantly a mixture 
of He, H, and with f >. 0.1, T < 100°K, and where cosmic rays make some 
contribution to the ionization rate. 

In Table I we list the main reactions which determine the ion and HD 
abundances in the clouds under consideration. The rates for cosmic ray 
ionization of H, and He include direct ionization by the cosmic rays 
and ionization by secondary electrons (Glassgold and Langer, 1973b, 1974). 

4 * 4 * 

The charge exchange reactions remove and return the ions H , D , and 0 . 

+ + 

Next listed are the ion-molecule reactions which make , HD, and OH ; 

OH subsequently reacts to fprm other oxygen-bearing molecules and molecular 
ions. The dissociative charge exchange reaction of He + with H^ dominates 
over other H -producing reactions with molecules because the large 
abundance of H^ more than compensates for the small rate constant k^. . 

Thus the fast reaction C + + OH -► CO + H + (Herbst and Klemperer, 1973; 

O’Donnell and Watson, 1974) cannot compete with k^ until the fractional 
abundance of OH becomes as large as 10 ^ £(He + )/E(C + ), where £(He + ) and 

4 * 4 * 4 * 

£(C ) are the fractional abundances of He and C . For diffuse and 
moderately thick clouds, this would require £(0H) 10 , which is 
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+ 

unreasonably large. We have also omitted other H -destroying charge exchange 
reactions with neutral molecules for similar reasons, i.e. because the rele- 
vant molecular abundances are too small in the clouds under consideration. 

For example, the rate constants for H + + OH H + 0H + and H + + H^O -*■ H + 

-93-1 + 

are probably 'v 10 cm s , and thus these reactions destroy H at a rate 

10-100 times slower than does charge exchange with 0 for temperatures in the 

range from 50-100°K. Our estimates on the abundances of the various oxygen 

molecules are based on a related investigation which we are now preparing 

for publication. In this work we also find that the densities of molecular 

ions such as OH , H^O , and H^O are negligible compared to the total ion 

density. Charge exchange and other ion-molecule reactions initiated by other 

. + + 

atomic ions such as C , He , etc. can be ignored for similar reasons. 

^ ^ -f. 

The chemical balance equations for 0 , D , He , , and H yield 

the following relations: 


n(Ee + ) = — SL ge ) S(He)n 

a(He )n g + H 2 ) 


n(0 + ) 


k 5 S(0)n n(H + ) 
k 6 n(H) + kgn(H 2 ) 


( 6 ) 


( 7 ) 


+ k 3 ?(D)nn(H + ) 

n ^ D 1, ’ k^n(H) + k^nCH.,) 


( 8 ) 



r 


n(H 3 + ) 


k 2 C 2 (H 2 + ) n 2 (H 2 ) 


e( H 3 + ) HgC^nCH) + k 2 n(H 2 )J 


and, for E , 


n(H + ) * - V(He + ) n(H 2 ) 

oc(H ) n + v n 
e 


( 10 ) 


>/here 




*l f / 2 


V^V 2 -*i> J 


(1 - f) + Co(H + )f/2, (11) 


8 


2 k + (k“ - 2k, ) ~ " k 5 ? <°> f ' 


( 12 ) 


and C(He), 5(0), and 5(D) are total fractional abundances relative to 
hydrogen. (Thus we have anticipated the conclusions implied by (6) - (8) 
that He, 0, and D are mainly neutral for the clouds under consideration. ) 

On the other hand, we take the elements with I < 13.6 eV to' be predominantly 
ionized, which requires N < 7 x 10 21 cm‘ 2 , assuming the ionizing radiation 
to be attenuated by a grain cross section of 1.25 x 10 -21 cm 2 per interstella: 

hydrogen atom. If 5 i is the total relative abundance of these elements, the 
total heavy Ion density is 

n i * ^i n ' 

and the electron density is 


: /‘T * L’J "-Xt v -y y f < ■ p 
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n = n(E + ) + n(He + ) + n. 
6 i 


( 13 ) 


The relative abundances of and , as well as all other atomic and 

—6 

molecular ions, are too small ( ^ 5 x 10 ) to be included in ( 13 ). 

An equation for the fractional abundance of electrons, x * n /n, is 

e e 

obtained by substituting (6) and (10) into (# ), 


e i 


i-)x - ( 

a e 


X 

a 


V 


V . + — 2 - 
i ii a 


■) = 0 , 


(14) 


where a * a(H ) and 


V ± - £. + x(He + ) 
and 

?' H “ Sjj + k x(He + )n 


(15) 


(16) 


Equation (14) is a cubic because of the a(He + ) x term in the denominator 

e 

of (6), but it can be simplified in a number of physically interesting 
cases. 

When ot(He ) f/2 , (14) becomes a quadratic with the 

solution: 

x e - 1/2 [(5'i - -£W (5'j + "a 1 ") 2 + 4 C' H /ha J . (I.) 
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x(H + ) - 1/2 [- (5' ± + -£-) + /(5'j + -*-) 2 + 4C' H /ncx J . (18) 

Using the values for a(He ) and h in Table I, the condition for ignoring 

7 

radiative recombination of He + is f » 200 x^ for temperatures in the range 

from 50-100°K. Since x is roughly of the order of 10 ^ provided 
^16 —X 

Cp ^ 10 s , (17) and (18) will apply as soon as a moderate fraction 

of hydrogen has been converted to molecular • form , say f = 0.1. 

When He radiative recombination can not be ignored, the ionization 
equations can be solved by iteration on the following system of 
equations: ‘ - 

z e “ [j^i “ + + + 4U f /na + x~ (He*)' ) J (19) 

,<»> 

e 9 

‘ ■ *. * i 

x(H + ) « x - x(He + ) - C, , . . (21) 

where 

j ‘ * 

?' « Cjj + ?(He ) c(He) (22) 

and where we have used a(He + ) * a(H + ) * a. Eq. (18) suggests another 

* r 

approximate solution, 



(23) 


x e 4 1/2 [ (? i - -|- )+ / <«,, + -J->* + 4 C'/na J 

i(H ) £ 1/2 (5^ + ) +/ (5^ + -J -) 2 + 4 C/na J (24) 

appropriate for x(He ) « l'/ny. Referring to (22) and (11) for ?' and 

(12) for y> we can see that this approximation tends to improve as T and 

£(0) are decreased. In fact this approximation works very well over the 

ranges of parameters considered in this work and is generally superior to (17) 
and (IS). • 

Electron and ion abundances are plotted in Figure 1 as a function of 

f for a typical parameter set: n » 100 cmf^, T = 8Q°K, £ * lCf^s - ^, 

£(0) = 10 , ^ * 7 x 10 (corresponding to heavy element depletion 

by a factor 7). The He fractional abundance decreases rapidly with " 
f (increasing abundance) as & result of dissociative charge exchange. 
Because this reaction produces H + ions, x(h + ) does not decrease rapidly 
until f becomes large (— 0.8). phe rapid decrease Occurs because of (i) the 
disappearance of direct cosmic ray ionization of H (c.f. the first term in 
(11) ), (ii) the smallness of dissociative ionization of (the second 

term in (ll) ), and (iii) 0 + ** OH + H competes with 0 + + H 0 + H + . 

The changes induced in Figure 1 by changes in parameters can be understood 
from (23) and (24) in terms of the two most important variables, £ /na 
and y * 10“ 9 e _232/T 5(0 )f cm^s** 1 . 

The appearance of in Figure 1 for large f is important because 
H^ + enters into numerous ion-molecule reactions (Herbst and Klemperer, 1973; 
Watson, 1973). In the context of the present discussion, H + reacts 
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exothermically with 0 to produce 0H + and at a rate % lO^cm^s” 1 , and 

ihus competes with endothermic H +0 charge exchange (k^ 'v 

when f £ 0.8 and T £ 100°K, even though it is only ^ 1/10 as abundant as H + . 

The relative importance of H + and E^ + depends on n, T, and £(0). For 

example, in a denser region where n = 10^ cm"^, T * 80°K, C * lO^^s -1 , 

. P 

and 5(0) « 10 , x(H^ ) * 0.1 x (H ) for f * 0.7 and x(H^ + ) * x(H + ) for 

x * 0.95. At lower temperatures, will be even more important. 

The reduction in the He abundance by dissociative charge exchange 
with has other implications for the ion-molecule reaction chemistry of 
the interstellar gas. Herbst and Klemperer (1973) assumed that He + is 
destroyed in dense clouds by dissociative charge excha ng e with CO, 

He + CO -*■ C + + 0 + He, (25) 

with a rate constant 2 x 10 7 cnTs (Laudenschlager and Bowers, 1973). In 
order for (25) to compete with (l), the relative CO abundance must be ^ 10~^. 
For smaller CO concentrations, the rate at which He + dissociatively charge 
exchanges with heavier neutral molecules, and thus produces atomic ions, 
is smaller than previously assumed. 

i 

III. Discussion. 

The abundance of HD provides a good example of the application of 
the results of Section II. Unless f is very' small (£ 10 “^), reaction (5) 
is much more important than grain formation for producing HD. The primary 
destruction mechanism is photo dissociation 
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■h'v + HD H + D , 


(26) 


> 


+ 

with rate constant n* The balance equations for HD and D (Black and 
Dalgamo, 1973; Jura, 1974; O’Donnell and Watson, 1974) lead to the 
relation 


x(HD) ‘ ~2 k, V - 2 k 4 )f *»>' f 


(27) 


We use the results of Section II to obtain the curves of Figure 2, using 
the parameters n * 10 cm , £(D) * 1.5 x 10 (Rogerson and York, 1973 )> 

£(0) * 10 ^ y “ 7 2 10 ^ , and n * 10 This value for n is based 

on the calculations of Stephens and Dalgamo (1972) and the average 

interstellar radiation field given by Jura (1974)- The solid curves 

o ~17 —17 

are for SO It and a range of primary ionisation rates, £ = 10 , 5 x 10 , 

—16 —1 

and 10 s , The existence of a strong maximum for x(HD) at an intermediate 

value of f occurs from the competition between the f- increasing coefficients 

in (27) and x(H + ), which Is a strongly decreasing function of f as f + 1. 

In applying these curves to cloud models, variations in physical properties 

of the clouds must be taken into account, such as density, temperature, 

radiation field, and size. For example, grain attenuation of the UV 

o 

radiation field near 1000 A reduces the dissociation rate n, roughly as 

?n —9 — 1 

n * n o e t = K(8 1 10 cm ) . (28) 


; + 

This decrease in T) tends to compensate for the decrease in x(H ), and 
thus to extend the region in which HD occurs with substantial abundance. 
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The two dashed curves in Figure 2 indicate the effects of varying 

— 1 7 —1 

the temperature for the particular case g * 5 x 10 s The values 

p 

60, 80, and 100°K roughly cover the range in kinetic temperatures observed 
for diffuse clouds (Spitzer and Cochran, 1973; Spitzer et al, 1974). The 
HD fraction is seen to change by as much as a factor of two for large f . 

The main source of the temperature variation is the e -2 ^ 2 ^ factor in 
the H + + 0 charge exchange rate constant k^ occurring in y (Eq. (12) ), 
which appears in the x(H + ) formula (24). 

In order to discuss the HD column density measurements of the 
Princeton group (Spitzer et al, 1973 and 1974), it is useful to rewrite 
(27) in terms of the HD-H 2 fraction y 5 n(HD)/n(H 2 ), 

y = ' 2k 4 "+ - & 4 )f k 3 5(13) ~ (p h ) ” * (29) ^ 

The advantage of (29) is that most of the variations in Ifche members of the 

right side tend to cancel for diffuse clouds. This can he seen in Figure 3, 

where y is plotted as a function of f for a number of chmces of c , 

P 

n = 20 and 100 cm, and for common values T = 80°K, £(G) = Kf'S 
-7 -10 -1 

= 5 x 10 , and n “ 10 s . We note that y is roughly constant 

-19 3 -1 

for f < 0.6 and C^/n > 5 x 10 cm s . Furthermore, when y is changed 
from 20 to 100 cm (corresponding to a reasonable density range for 
diffuse clouds, Glassgold and Langer, 1973a, 1974; Jura, 1974), y changes 
by at most a factor of 2. According to our earlier model calculations, n 
does not vary by more than a factor of 2 and T changes little in going 
from the outside to the inside of a diffuse cloud. Thus we can ignore 
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the variation of y with n and T. On the other hand, UV attenuation (28) 

can be important, depending on the column density N into the cloud. The 

effect of decreasing T) may be crudely viewed as a raising of the large f 

portions of the curves of Figure 3 so as to maintain y roughly constant 

up to f 0.8 - 0.9. To summarize, (29) together with the ionization 

formulae (24), provide justification for replacing y « rt(HD)/n(H 2 ) by 

the column density ratio Y = N £kd] /N > as was done in previous work 

(Black and Dalgamo, 1973; Jura, 1974; O'Donnell and Watson, 1974). The 

measurements of Y (Spitzer et al, 1973, 1974) vs. F * 2M [h 2 3 /N have 

also been included in Figure 3 mainly for reference, and not for purposes 

of direct comparison with the theoretical curves of y vs. f, A proper 

comparison can only be made in terms of a more complete (inhomogeneous) 

model calculation for Y vs. N, such as we recently reported for H 2 and 

other species (Glassgold and Langer, 1974). On the other hand, it is 

interesting that a number of the measurements of Y cluster about the value 
-7 

8 x 10 , and that the deviations for 59 Cyg and % Qph may be partly 

understood in terms of variations in the relevant radiation field (Jura, 
1974). O'Donnell and Watson (1974) have proposed a similar explanation 
for the anomalously large value of Y measured for 10 Lac. 

Proceeding on the basis that y « Y ■ constant for a diffuse cloud, 
we can invert (29) and obtain a rough formula for the dependence of n(H + ) 
on N and f , 

n(H + ) = £y\ + (V2 - VV f V ~ N(8 X 1 ° 2 ° cm ’ 2) ’ 1 Y. (30) 
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According to this result, substantial reductions are possible for n(H + ) in 
going from the outside to the inside of a cloud, in part from attenuation 
and in part from the coefficient in square brackets. The numerical values 
for n(H ) obtained from (30) at the center of a diffuse cloud are essentially 
the same as those obtained by previous authors, i.e. in the range from 10"^ - 
lCf 2 cm -3 . 

To conclude this discussion of the HD problem, it must be said that 
a quantitative understanding of the observations must still be given. We 
have attempted to emphasize the importance of various inhomogeneities, 
especially the variation in i(E + ) and the attenuation of the HD photo- 
dissociation rate. We did find some basis for the possibility that these 
two effects partly compensate in diffuse clouds, but a more detailed model 
calculation is really needed. The results of any calculation will be 
strongly influenced by the radiation field appropriate to each cloud and 
by the values for rather uncertain charge- exchange and ion-molecule 
reaction rate constants. The fact that some of the measured values for 
Y are much larger or much smaller than the average may prove difficult to 
explain. In view of all these qualifications and uncertainties, we believe 

that the present HD abundance determinations cannot determine £ to much 

P 

better than a factor of 5. As previous authors have estimated, r is 

P 

likely to lie somewhere in the range from 1 CT's~' 1 ' - lCT^s"^. 

The results of Section XI also relate to the determination of inter- 
stellar cloud electron densities from absorption line measurements. Since 
we have recently commented on observations of Cal and Call (White, 1973; 
Bortolot et al, 1974) and Cl and CII (Morton et al, 1973), we note that Hobbs 
(1973, 1974) and Lutz (1974) have recently made deductions about ionization 
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in clouds from KI absorption lines. Hobbs attempts to draw conclusions 

about the density variation of x and the difference between x and x(H + ). 

e e 

-I 

He finds that the n variation of x associated with cosmic ray ionization 

e ^ 

of atomic H is in conflict with KI measurements, and he concludes that 
—16 —1 

Sp a. 10 s The theory of Section II contains a number of improvements 

in previous discussions of ionization in clouds which are relevant in this 

connection. In particular, we have included the important H + destruction 

mechanism of charge exchange with 0, the suppression of direct ionization 

of H by the formation of molecular hydrogen, and the transformation of He + 

into H ions by dissociative charge exchange. The latter effects are 

contained in the effective ionization rates and £ r defined by (ll) 

and (22), respectively. Taken together with the above discussion of HD 

(which suggests that c lies in the range lCf 17 - 1(T 16 s" 1 ), it does seem 

* 

on the basis of Eq, (24) that a clear n ^ variation for x will not be 

e 

manifested in interstellar clouds except for very small densities. 

In conclusion, we find that the ionization structure of a diffuse 
cloud containing moderate amounts of molecular hydrogen is inhomogeneous. 
Even for the clouds observed by the 0A0-3 program, fairly large values of 
f are appropriate toward their centers, so that much of the variation 
in d i cated in Figure 1 may be expected to occur. Additional inhomogeneities 
arise from the density and tenperature variation predicted by model 
calculations (Glassgold and Langer, 1974). 

This work has been supported by Grant NGR- 3 3-016-196 from the 
National Aeronautics and Space Administration. • 
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Table I. 

.Reaction 

' . Rate 

Reference 

Cosmic Ray Ionization 

. 

1 

„ + 


' 7 

p + H H + e + p 



p + + e + p 

? 2^ H 2 ^ = 2,3 


p + H 2 **H + H + + e+ p 

C 2 (H + ) » 0.02 ? 2 (H 2 + ) 


+ 

p + He -*■ He + e + p 

?(He) = 1.5 C 


Charge Exchange 

' 


•;>, h 2 + + h -j- h 2 + H + 

kj- = 6 x 10" 10 cm 3 s -1 

2 

* D + H + -»■ D + + H 

• k 3 = 2 x IQ" 9 e-^cmV 1 ; 

3 

D + + H D + H + 

k, = 2 x lO^omV 1 

3 

0 + H + -*• 0 + + H 

: k = lo-V 232 ^- 1 ' 

4 

0 + + H + 0 + H + 

^ = lO^cmV 1 

; 4 

Ion Molecule Reactions 



: r 2 * + h 2 + H + + H 

k 2 * 2 x 10“ 9 cm 3 s' 1 

5 

’ D + + Hg -*• HD + H + 

ky = 10 _9 cm 3 s _1 

6 

0 + + H 2 0H + + H 

kg = 2 x 10 -9 cm 3 s -1 

7 

Dissociative Charge Exchange 


■ 

He + ** He + H + H + 

kg = 10“ 13 cm 3 s -1 

8 

• k 

Radiative and Dissociative Recombination >- 




(Table 1, contd. ) 


Reaction 

+ 

e + He -*• He + y 
e + H 2 + 3H or H 2 + H 


!• Glassgold and Langer (1973) 

2. de Jong (1972) 

3. Smith (1966) 


4* Field and Steigman (1971) 

5. Neynaber and Trujillo (1972) 

6. Fehsenfeld et al (1973) 

7. Fehsenfeld et al (1966) 

S. Biondi (1974) 

9. Spitzer (1966) 

10. Leu, Johnson, and Biondi (1973) 



FIGURE CAPTIONS 


Fractional abundances of electrons and various ions as a function of 
fractional abundance f * 2n(H 2 )/n, n - n(H) + 2n(H 2 ) for T “ 80°K, 

n - 100 cm” 3 , c = 10“ 16 s _1 , £(0) = 10"*, and £ . = 7 x 10~ 5 . The 

Jr ^ 

dashed line is the fractional abundance of ions of atoms 

with ionization potential I < 13.6 eV. 


Fractional abundance of K + ions as a function of f for T = 80°K, 

n = 100 cm -3 , 5(0) = 10 -4 , q ■ 7i 10 -5 and r = 10' 17 , 5 x 10 -17 , 

“16 “1 o 

and 10 s (solid curves). The dashed curves are for T = 60 U K and 

100°K and C =5x 10“ 17 s _1 . 

P 


Ratio of ED to densities, y = n(HD)/n(H 2 ) as a function of f for 

T ~ 80°K, £(0) = 10^, £. = 7 x 10 and a range of primary cosmic 

ray ionisation rates t, . The solid curves are for a density n = 100 cm , 

P 

-3 ' 

and the dashed curves for n = 20 cm . The open circles are the measured 
ratios of Y = N [wJ/N as a function F * 2N /N, N = n£hJ+ 2R C\~2 

(on the same • scale as y and f, respectively) as reported by Spitzer 


et al (1973, 1974). 
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